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PRELIMINARY AMENDMENT 

Assistant Commissioner for Patents 
Washington, D.C. 20231 

Sir: 

Prior to examination of the above-captioned patent application, applicant requests 
that the following claim amendments be entered. Additional amendments are incorporated 
in the formatted Substitute Specification submitted herewith. A copy of the English- 
language translation of the original application, together with a CompareRite® version 
showing the amendments made to the original application, in accordance with 37 C.F.R. 
§1. 121 (2001), are also submitted herewith. No new matter has been introduced in these 
amendments to the original specification. 

IN THE CLAIMS: 

Please replace Claims 1-9, as follows: 

1. (Amended) A method for optimizing power plant control values, where the 
method comprises the steps of 
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a) determining a trajectory of plant control values, 

b) determining from the trajectory of plant control values, by simulation of the 
power plant's behavior, corresponding trajectories of costs for operating the 
plant and trajectories of revenues from selling produced output, 

c) computing a total operating cost as an integral of an objective function that 
comprises a difference between costs and revenues, 

d) iteratively repeating steps a) through c) with the optimization module varying 
the trajectory of plant control values until an optimized lower total operating 
cost is arrived at, 

wherein in step b), by modeling of plant components, a trajectory representing a rate of 
ageing of plant components is determined from the trajectory of plant control values, and 
that, in step c), the objective function comprises said rate of ageing. 

2. (Amended) Method according to claim 1, wherein, in step b) trajectories of 
plant output values are determined, and that, in step c) the objective function comprises a 
difference between plant output values and associated demand values. 

3. (Amended) Method according to claim 1, wherein the total operating cost is 
computed as 
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where R l5 R 3 are weighting matrices, t is a present time and T is a predetermined duration. 

4. (Amended) Method according to claim 2, wherein the total operating cost is 
computed as 

t + T 

Jlu] = f R x e(x) + R 2 [d(x)-p(x)] 2 + R 3 [c(x)-q(x)] dx, 

where R 19 R 2 , R 3 are weighting matrices, t is a. present time and T is a predetermined 
duration. 

5 . (Amended) Method according to claim 1 , wherein at least part of the 
trajectory of plant control values corresponding to the optimized lower total operating cost 
is displayed to an operator. 

6. (Amended) Method according to claim 1, wherein at least part of the 
trajectory of plant control values corresponding to the optimized lower total operating cost 
is used to control the actual plant 3. 

7. (Amended) A computer program product comprising a computer readable 
medium, having thereon: computer program code means to make, when said program is 
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loaded in a computer that is operationally connected to a power plant, the computer execute 
the method according to claim 1 . 

8. (Amended) A power plant optimizing system that optimizes plant control 
values and comprises 

a) an optimization module that is configured to determine a plurality of 
trajectories of plant control values and to determine a total operating cost as an 
integral of an objective function that comprises a difference between costs and 
revenues that are associated with operating the power plant according to the 
plant control values, 

b) a simulation module that is configured to determine, from the trajectory of 
plant control values, by simulation of the power plant's behavior, 
corresponding trajectories of costs for operating the plant and trajectories of 
revenues from selling produced output, 

wherein the simulation module is configured to determine, from the trajectory of plant 
control values, a trajectory representing a rate of ageing of plant components, and that the 
objective function comprises said rate of ageing. 
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9. (Amended) A power plant optimizing system according to claim 8, wherein 
the simulation module is configured to determine, from the trajectory of plant control 
values, plant output values, and that the objective function comprises a difference between 
plant output values and associated demand values. 
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REMARKS 



By way of the foregoing amendments to the claims, Claims 1-9 have been amended 
to delete the multiple dependencies and reference numerals, and to replace the words 
"characterized in that" with the word "wherein". These changes have been made in 
accordance with 37 C.F.R. § 1. 121 as amended on November 7, 2000. Marked-up 
versions of Claims 1-9 indicating the changes accompany this Preliminary Amendment. 

Early and favorable consideration with respect to this application is respectfully 
requested. 

Should any questions arise in connection with this application, the undersigned 
respectfully requests that he be contacted at the number indicated below. 



Respectfully submitted, 

Burns, Doane, Swecker & Mathis, l.l.p. 



By: 




Robert S . Swecker 
Registration No. 19,885 



P. O. Box 1404 

Alexandria, Virginia 22313-1404 
(703) 836-6620 



Date: December 19, 2001 
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Attachment to Preliminary Am endment dated December 19. 2001 
Marked-up Claims 1-9 

1. (Amended) A method for optimizing power plant control values, where the 
method comprises the steps of 

a) determining a trajectory of plant control values [u(x)], 

b) determining from the trajectory of plant control values [ u(x) ] , by simulation of the 
power plant's behavior, corresponding trajectories of costs [c(x)] for operating the 
plant and trajectories of revenues [q(x)] from selling produced output, 

c) computing a total operating cost [JfuJ] as an integral of an objective function that 
comprises a difference between costs [c(x)] and revenues [q(x)], 

d) iteratively repeating steps a) through c) with the optimization module varying the 
trajectory of plant control values [ u(x) ] until an optimized lower total operating 
cost [JfuJ] is arrived at, 

[characterized in that] wherein in step b), by modeling of plant components, a trajectory 
representing a rate of ageing [e(x)] of plant components is determined from the trajectory of 
plant control values [u(x)], and that, in step c), the objective function comprises said rate of 
ageing [e(x)]. 

2. (Amended) Method according to claim 1, [characterized in that] wherein , in step 
b) trajectories of plant output values [p(x)] are determined, and that, in step c) the objective 
function comprises a difference between plant output values [p(x)] and associated demand 
values [d(x)]. 
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Attachment to Preliminary Ampnriment rlat e d December 19. 2001 
Marked-up Claims 1-9 

3. (Amended) Method according to claim 1 , [characterized in that] wherein the total 
operating cost [JfuJ] is computed as 

ti-T 

J[u] = f R x e(x) + 7?3[c(t)-^(x)] dx, 

t 

where R 1? R 3 are weighting matrices, t is a present time and T is a predetermined duration. 

4. (Amended) Method according to claim 2, [characterized in that] wherein the total 
operating cost [JfuJ] is computed as 

t+T 

J[u] = J R,e(x) + R 2 [d(x)-p(x)] 2 + R 3 [c(x)-q(x)] dx, 

where R v R 2 , Rj are weighting matrices, t is a present time and T is a predetermined duration. 

5. (Amended) Method according to claim 1, [characterized in that] wherein at least 
part of the trajectory of plant control values [u(x)~\ corresponding to the optimized lower total 
operating cost [JfuJ] is displayed to an operator. 
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Attachment to Preliminary Am endment dated December 19. 2001 
Marked-up Claims 1-9 

6. (Amended) Method according to claim 1, [characterized in that] wherein at least 
part of the trajectory of plant control values [w(x)] corresponding to the optimized lower total 

m operating cost [J[uJ] is used to control the actual plant 3. 

7. (Amended) A computer program product comprising a computer readable 
j: medium, having thereon: computer program code means to make, when said program is loaded 
H in a computer that is operationally connected to a power plant, the computer execute the 

ry 

=T method according to [one of claims 1 to 6] claim 1 . 

8. (Amended) A power plant optimizing system that optimizes plant control values 
and comprises 

a) an optimization module [(1)] that is configured to determine a plurality of 
trajectories of plant control values [u(x)] and to determine a total operating cost 
JfuJ as an integral of an objective function that comprises a difference between 
costs [c(0] and revenues [q(x)] that are associated with operating the power plant 
according to the plant control values [ u(x) ] , 

b) a simulation module [(2)] that is configured to determine, from the trajectory of 
plant control values [u(x)], by simulation of the power plant's behavior, 
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Attachment, to Pr eliminary Amendm ent dated December 19. 2001 
Marked-up Claims 1-9 

corresponding trajectories of costs [c(x)] for operating the plant and trajectories 

of revenues [ q(x) ] from selling produced output, 
[characterized in that] wherein the simulation module [(2)] is configured to determine, from 
the trajectory of plant control values [u(x)], a trajectory representing a rate of ageing [e(x)J 
of plant components, and that the objective function comprises said rate of ageing [e(x)]. 

9. (Amended) A power plant optimizing system according to claim 8, [characterized 
in that] wherein the simulation module [(2)] is configured to determine, from the trajectory of 
plant control values [w(t)], plant output values [p(x)], and that the objective function 
comprises a difference between plant output values [p(x)] and associated demand values 
[d(x)]. 



COMPARERITE® SPECIFICATION 
(WITHOUT CLAIMS) 



j-Field of the Invention 

[[0001]] The invention relates to the field of process control, in particular to a method and 
computer program product for optimizing power plant control values and a power plant 
optimizing system as described in the preamble of claim 1, 7 and 8, respectively. 

Background of the Invention 

[[0002]] Methods and systems for optimizing power plant operation are shown in US 
5,347,466 and US 5,886,895. In both cases, plant operation, which is controlled by 
control values generated by a control system, is optimized in order to minimize fuel costs 
while achieving predetermined required output values for produced power and process 
steam. This is done by, at a given time, determining future values of control values and 
simulating a behavior of the plant up to a given future time. From the simulation, fuel 
costs are determined and a objective function comprising the fuel costs and costs for buying 
power from another source is computed. In order to determine optimal control values, the 
control values are varied and the simulation is repeated until a minimum of the objective 
function is found. However, such an optimization is constrained by the fact that 



predetermined operating limits of the plant may not be exceeded. This has the effect that 
the optimization does not necessarily give a financially optimal result for the plant owner. 



{Description} 

rSiimmaryl of the Invention 

[[0003]] It is an object of the invention to create a method and computer program product 
for optimizing power plant control values and a power plant optimizing system of the type 
mentioned initially that overcomes the deficiencies mentioned above. 
[[0004]] These objects are achieved by a method and computer program product for 
optimizing power plant control values and a power plant optimizing system according to the 
claims 1, 7 and 8. 

[[0005]] In the inventive method and computer program product for optimizing power 
plant control values and a power plant optimizing system, a rate of ageing of plant 
components is determined when simulating a future behavior of the plant, and an objective 
function to be minimized by the optimization comprises said rate of ageing. 
[[0006]] In this way, the invention incorporates long term financial effects of component 
ageing into short term decisions on plant operation. For example, under certain conditions 
and when energy prices are high, the optimization shall indicate that it is financially 
advantageous to operate the plant outside normal conditions, i.e. that higher revenue from 
generated power more than compensates for higher maintenance or replacement costs due 
to increased component degradation. 



[[0007]] In a preferred variant of the invention, plant output values representing plant 
production are determined in the simulation and a difference between said plant output 
values and demand values is included in the objective function. 

[[0008]] Further preferred embodiments are evident from the dependent patent claims. 
Brief Description of th e Drawings 

[[0009]] The subject matter of the invention will be explained in more detail in the 
following text with reference to preferred exemplary embodiments which are illustrated in 
the attached drawings, in which: 

[[0010]] Figure 1 shows a block diagram corresponding to the invention; and 
[[0011]] Figure 2 shows a block diagram corresponding to a preferred embodiment of 
the invention. 

[[0012]] The reference symbols used in the drawings, and their meanings, are listed in 
summary form in the list of reference symbols. In principle, identical parts are provided 
with the same reference symbols in the figures. 
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Detailed D e scription of (Preferred Embodimen ts 
}[the Invention 

[0013]] Figure 1 shows a block diagram corresponding to the invention. An optimization 
module 1 is provided with demand valuesfr^- d(x). The optimization module 1 generates 
power plant operating parameters or control values u(x) that are input to a simulation 
module 2. The simulation module 2 comprises a simulation model of an actual power plant 
3 and determines values for a power plant output or productionfr} p(x), revenue values 
q(x) from selling the produced output, costs c(x) associated with a given power plant 
output, and values of a rate of ageing e(x) of plant components. The values determined by 
the simulation module 2 are input to the optimization module 1 . 
[[0014]] The values d(x),u(x),p(x),^(x),c(x),e(x) mentioned above are, in general, 
vectors that change over time, i.e. each vector represents a trajectory in a multidimensional 
space. For example, d(x) represents an changing amount of electric power and of process 
steam to be generated over the next hour or day. 

[[0015]] The control values u(x) represent plant input or operation parameters that are used 
to control the plant, such as valve and guide vane positions, or set points for local 
controllers of temperature, fuel flow or steam flow. From a trajectory of control values 
u(x) that is given for a time interval from a present time [t] to a future time ft + 7], where 
[T] is called "prediction horizon" or "optimization horizon", the simulation model 
determines a trajectory of the state of the power plant, i.e. the values of a plurality of 
variables that represent the state of the plant as it changes with time. From the state 
trajectory, trajectories p(x),q(x), c(x),e(x) for the values output by the simulation module 2 
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are determined for said time interval. In a preferred embodiment of the invention, the 
simulation model is a dynamic model of the power plant and its components. 
[[0016]] The production p(x) represents the amount of electric power and process steam 
that would be generated if the control values «(T)were applied to the actual power plant 3. 
The revenue term q(x) represents revenues from selling the electric power and process 
steam and is computed by multiplying the production with a price per unit of production. 
Computing the revenue q\x) preferably includes power prices that vary over time, either 
because agreements for providing a given amount of power at a given price and at a given 
time exist, or because estimates about future spot market power prices are available. 
[[0017]] The costs c(t) are computed from the fuel useage, fuel costs, and process 
efficiencies required to produce, p(x) . The fuel consumption and process efficiencies are 
computed using a plant simulation model and depend directly on the plant control values 
it(x) . Associating fuel costs with outputs is done, for example, by 

m 

[[0018]] assigning all fuel costs to electric power output and considering the produced 
steam as a free by-product, and 

m 

[[0019]] assigning additional fuel costs to the steam output, if they arise because a higher 
steam production is required than would be produced at the required level of electrical 
power 

[[0020]] The rate of ageing e(x) represents costs associated with a reduction of lifetime of 
plant components. Plant components are e.g. pressurized vessels, pipes or turbine blades. 
They are designed to withstand normal operating conditions for a given amount of time, 
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i.e. their lifetime under nominal conditions. Subjecting a component to more severe 
operating conditions such as higher temperatures or temperature gradients increases 
component degradation and reduces component lifetime. This in turn increases long term 
costs for operating the power plant. According to the invention, an effect on these long 
term costs is included in the total operating costs. 

[[0021]] The rate of ageing e(i) is computed from models associated with various plant 
components, such as the following models, which are well know to one versed in the art: 

m 

[[0022]] Microcrack growth rate models for components operating at elevated 
temperatures, such as boiler pressure parts, steam pipes and steam turbine rotors. See for 
example "Life Assessment Technology for Fossil Power Plants", R. Viswanathan, Proc. 
American Power Conf., Vol. 57-111: 1749-1768, 1995 and "A Lifetime Prediction Model 
for Thermomechanical Fatigue based on Microcrack Propagation" , {M.P.Mill e r et.al} 
[M.P. Miller et. al], ASTM STP 1186, 1993, pp. 35-49. 

m 

[[0023]] Corrosion, oxidation, and erosion models based on extensions to microcrack 
growth models as described in "Corrosion Fatigue: Mechanics, Metallurgy, 
Elctrochemistry, and Engineering", T. Crooker and B. Leis editors, ASTM 1983. 

m 

[[0024]] Mechanical and thermomechanic stress models to determine stresses created in 
plant components, e.g. stresses in a steam turbine rotor shaft due to load oscillations and/or 
steam temperature changes. 



[[0025]] The rate of ageing e(x) represents the plant lifetime consumption or an asset value 
degradation expressed as a financial value under operating conditions caused by the control 
values w(t) . The entries in e(z) correspond to the aging rate for each plant component and 
are, for example, computed directly from 

R da 



where [a crUJ is a predefined critical crack length for said component, [a] is the crack length 

computed using the models described above and dependent on u(x) , and, [R\ , is the 

component, cost, repair, or maintenance cost. 

[[0026]] In the inventive method, the following steps are performed 

[[0027]] a) The optimization module 1 determines the trajectory of plant control values 

u(z) from the present time ft] for a predetermined duration [T] and transmits said plant 

control values u(x) to the simulation module 2. The first time that this step is executed, 

plant control values are e.g. set to constant values for the duration [T] or are determined 

from stored standard trajectories corresponding to a given operating condition such as a 

plant startup or a load change. Later executions of this step preferably use a variation of 

the plant control values u(x) of an earlier step. 

[[0028]] b) The plant model incorporated in the simulation module 2 determines from 
the trajectory of control values u(x) , by simulation of the power plant behavior, a 
trajectory of the state of the power plant. From said state, corresponding trajectories of 
costs c(x) for operating the plant and trajectories of revenues q(x) from selling produced 



output are determined, as well as the rate of ageing e(x) of plant components. In a 
preferred embodiment of the invention, trajectories of plant output values p(x) are 
determined as well. 

[[0029]] c) A total operating cost [J[uJ] is computed as an integral of an objective 
function that comprises the rate of ageing e(x) and a difference between costs c{x) and 
revenues q(x) . In a preferred embodiment of the invention, this total operating cost [JfuJJ 
is defined as the functional 



where [R 1? R 3] are weighting matrices. The weighting matrices determine a relative 
influence of the magnitude of each component of each of the different vectors 
q(x),c(x),e(x). 

[[0030]] In a preferred embodiment of the invention, the above objective function also 
includes a difference between plant output values ^(-c) and demand values d(x) . In this 
embodiment, the total operating cost [JfuJJ is preferably defined asfr}- 



where [R t , R 2 , R 3] are weighting matrices. The terms in the objective function are 
interpreted as follows. The u [d(v) - p(v)Y term represents a penalty associated with not 




J[u] = f R.eix) + R 2 [d(x) -p(x)f +i? 3 [c(x)-q(x)] dx, 



satisfying the required power output. The "fc(r) - q(r)Y term represents costs minus 
revenue, i.e., an economic performance {indicator. d)} [indicator. 
[0031] d)] Steps a) through c) are repeated iteratively, with the optimization module 
varying the trajectory of plant control values [u:[t, t + T]»R nJ until an optimized lower 
total operating cost [JfuJJ is arrived at. As the functions q(x),q(x), c(x),e(x) depend upon 
plant control values u(x) , the optimization routine finds a strategy [u:[t, t + T]"R nJ that 
minimizes the total operating cost [JfuJ. 
[0032]]{- 

} Examples for units are, 

m 

[[0033]] for the rate of ageing e : $ per second 

m 

[[0034]] for costs c or revenue q : $ per kW or $ per kg steam per second 

m 

[[0035]] for Rl: dimensionless 

m 

[[0036]] for R2: $/kW A 2 or $(s/kg) A 2. 

m 

[[0037]] for R3: kW or kg per second 

m 

[[0038]] for demand d or production p : kW or kg per second 

[[0039]] In another preferred embodiment of the invention, the objective function is 
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ti-T 

J[u] = f R,m + R^d^-pix)? + R 3 c(x)-R 4 q(x) dz, 



where costs and revenues are weighted separately, using an additional weighting matrix 
[R 4 ], which allows for the cost and revenue vectors c(t) , q(x) to have different 
dimensions. 

[[0040]] The optimization step c) involves the problem of how to find the optimal solution 
efficiently, since the search space determined by the plurality of plant control values and 
the duration [TJ is large. Standard optimization procedures may be used, or dedicated 
solutions as described, for example, in US 5,347,466, US 5,886,895 or in "The 
Generalized Unit Commitment Problem", Ross Baldick, IEEE Transactions on Power 
Systems, 10(1): 465-475, February 1995. 

[[0041]] In a preferred embodiment of the invention, at time ft], the control strategy, i.e. 
the future values [u:[t, t + TJ»R n} corresponding to optimal plant operation, is displayed 
graphically to an operator, providing him with guidance for controlling the plant. 
[[0042]] Figure 2 shows a block diagram corresponding to another preferred embodiment 
of the invention, in which, at time ft], the control value [u(t)J is applied to the actual 
power plant 3, resulting in automatic optimized plant operation. Preferably, measured 
plant values y(t) of the plant are input to the simulation module 2 in order to update the 
simulation model. Appropriate modeling and model update techniques are well known 
from the field optimal control. In particular, receding horizon control techniques might 
result useful. The power plant is, for example, a gas turbine, a combined cycle power 
plant, or of a boiler-turbine-generator type. 
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[[0043]] A computer program product according to the invention comprises a computer 
readable medium, having thereon: computer program code means to make, when said 
program is loaded in a computer that is operationally connected to a power plant, the 
computer execute the method according to one of claims 1 to 6. A computer program 
according to the invention comprises said computer program code. 

[[0044]] A power plant optimizing system that optimizes plant control values according to 
the invention comprises 

[[0045]] a) an optimization module 1 that is configured to determine a plurality of 
trajectories of plant control values u(x) and to determine a total operating cost [JfuJJ as an 
integral of an objective function that comprises a difference between costs c(x) and 
revenues q(x) that are associated with operating the power plant according to the plant 
control values u(x) , and 

[[0046]] b) a simulation module 2 that is configured to determine, from the trajectory of 
plant control values u(x) , by simulation of the power plant's behavior, corresponding 
trajectories of costs c(x) for operating the plant and trajectories of revenues q(x) from 
selling produced output. 

[[0047]] The simulation module 2 is configured to determine, from the trajectory of plant 
control values ii(x) , a trajectory representing a rate of ageing e(x) of plant components. 
The objective function comprises said rate of ageing e(x) . 

[[0048]] In a preferred embodiment of the power plant optimizing system, the simulation 
module 2 is configured to determine, from the trajectory of plant control values u(x) , plant 
output values p(x) , and the objective function comprises a difference between plant output 
values p(x) and associated demand values d(x) . 
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[[0049]] The invention allows the operator or an automatic controller to incorporate not 
only short term financial considerations into a plant operation strategy, but also long term 
financial considerations that involve maintenance and replacement costs and their 
dependency on operating conditions. 
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List of designations 

[[0050]] 1 optimization {module!} [module 

simulation {modulc3} [module 
actual power plant 
control values 
production 
costs 
demand 
revenue 
rate of ageing 
measured plant values 
weight matrices 



[0051] 2] 
[0052] 3] 
[[0053]] u(x) 
[[0054]] p(x) 
[[0055]] c(x) 
[[0056]] d(x) 

[[0057]] m 

[[0058]] e(x) 
[[0059]] 
[[0060] R 1? R 2 , R 3 ,] 
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[ABSTRACT] 

In a method and computer program product for optimizing power plant control 
values and a power plant optimizing system an optimization module (1) minimizes total 
plant operation costs while achieving predetermined required output values for produced 
power and process steam. This is done by, at a given time, determining future values of 
control values and simulating, in a simulation module (2), the behavior of the plant up to a 
given future time. Corresponding fuel costs and generated power are determined in the 
simulation, and incorporated in an objective function. The optimization module (1) 
minimizes the objective function by varying the control values. According to the 
invention, a rate of ageing of plant components is determined when simulating the future 
behavior of the plant, and the objective function to be minimized comprises said rate of 
ageing. 
{(figure 1)} 
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Field of the Invention 

[0001] The invention relates to the field of process control, in particular to a method and 
computer program product for optimizing power plant control values and a power plant 
optimizing system as described in the preamble of claim 1, 7 and 8, respectively. 

Background of the Invention 

[0002] Methods and systems for optimizing power plant operation are shown in US 
5,347,466 and US 5,886,895. In both cases, plant operation, which is controlled by 
control values generated by a control system, is optimized in order to minimize fuel costs 
while achieving predetermined required output values for produced power and process 
steam. This is done by, at a given time, determining future values of control values and 
simulating a behavior of the plant up to a given future time. From the simulation, fuel 
costs are determined and a objective function comprising the fuel costs and costs for buying 
power from another source is computed. In order to determine optimal control values, the 
control values are varied and the simulation is repeated until a minimum of the objective 
function is found. However, such an optimization is constrained by the fact that 
predetermined operating limits of the plant may not be exceeded. This has the effect that 
the optimization does not necessarily give a financially optimal result for the plant owner. 

Summary of the Invention 

[0003] It is an object of the invention to create a method and computer program product 
for optimizing power plant control values and a power plant optimizing system of the type 
mentioned initially that overcomes the deficiencies mentioned above. 



[0004] These objects are achieved by a method and computer program product for 
optimizing power plant control values and a power plant optimizing system according to the 
claims 1, 7 and 8. 

[0005] In the inventive method and computer program product for optimizing power 
plant control values and a power plant optimizing system, a rate of ageing of plant 
components is determined when simulating a future behavior of the plant, and an objective 
function to be minimized by the optimization comprises said rate of ageing. 
[0006] In this way, the invention incorporates long term financial effects of component 
ageing into short term decisions on plant operation. For example, under certain conditions 
and when energy prices are high, the optimization shall indicate that it is financially 
advantageous to operate the plant outside normal conditions, i.e. that higher revenue from 
generated power more than compensates for higher maintenance or replacement costs due 
to increased component degradation. 

[0007] In a preferred variant of the invention, plant output values representing plant 
production are determined in the simulation and a difference between said plant output 
values and demand values is included in the objective function. 

[0008] Further preferred embodiments are evident from the dependent patent claims. 
Brief Description of the Drawings 

[0009] The subject matter of the invention will be explained in more detail in the 
following text with reference to preferred exemplary embodiments which are illustrated in 
the attached drawings, in which: 

[0010] Figure 1 shows a block diagram corresponding to the invention; and 
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[0011] Figure 2 shows a block diagram corresponding to a preferred embodiment of 
the invention. 

[0012] The reference symbols used in the drawings, and their meanings, are listed in 
summary form in the list of reference symbols. In principle, identical parts are provided 
with the same reference symbols in the figures. 

Detailed Description of the Invention 

[0013] Figure 1 shows a block diagram corresponding to the invention. An optimization 
module 1 is provided with demand values d(x). The optimization module 1 generates 
power plant operating parameters or control values w(x) that are input to a simulation 
module 2. The simulation module 2 comprises a simulation model of an actual power plant 
3 and determines values for a power plant output or production p(x), revenue values q(x) 
from selling the produced output, costs c(t) associated with a given power plant output, 
and values of a rate of ageing e(x) of plant components. The values determined by the 
simulation module 2 are input to the optimization module 1 . 

[0014] The values d(x), u(x),p(x), q{x),c(x), £(x) mentioned above are, in general, 
vectors that change over time, i.e. each vector represents a trajectory in a multidimensional 
space. For example, d{x) represents an changing amount of electric power and of process 
steam to be generated over the next hour or day. 

[0015] The control values u(x) represent plant input or operation parameters that are used 
to control the plant, such as valve and guide vane positions, or set points for local 
controllers of temperature, fuel flow or steam flow. From a trajectory of control values 
u(x) that is given for a time interval from a present time t to a future time t+T, where T is 



called "prediction horizon" or "optimization horizon" , the simulation model determines a 
trajectory of the state of the power plant, i.e. the values of a plurality of variables that 
represent the state of the plant as it changes with time. From the state trajectory, 
trajectories p(x), q(x), c(x), e(x) for the values output by the simulation module 2 are 
determined for said time interval. In a preferred embodiment of the invention, the 
simulation model is a dynamic model of the power plant and its components. 
[0016] The production p(x) represents the amount of electric power and process steam 
that would be generated if the control values w(T)were applied to the actual power plant 3. 
The revenue term q(x) represents revenues from selling the electric power and process 
steam and is computed by multiplying the production with a price per unit of production. 
Computing the revenue q(x) preferably includes power prices that vary over time, either 
because agreements for providing a given amount of power at a given price and at a given 
time exist, or because estimates about future spot market power prices are available. 
[0017] The costs c(t) are computed from the fuel useage, fuel costs, and process 
efficiencies required to produce, p(x) . The fuel consumption and process efficiencies are 
computed using a plant simulation model and depend directly on the plant control values 
u(x) . Associating fuel costs with outputs is done, for example, by 
[0018] assigning all fuel costs to electric power output and considering the produced 
steam as a free by-product, and 

[0019] assigning additional fuel costs to the steam output, if they arise because a higher 
steam production is required than would be produced at the required level of electrical 
power 



[0020] The rate of ageing e(x) represents costs associated with a reduction of lifetime of 
plant components. Plant components are e.g. pressurized vessels, pipes or turbine blades. 
They are designed to withstand normal operating conditions for a given amount of time, 
i.e. their lifetime under nominal conditions. Subjecting a component to more severe 
operating conditions such as higher temperatures or temperature gradients increases 
component degradation and reduces component lifetime. This in turn increases long term 
costs for operating the power plant. According to the invention, an effect on these long 
term costs is included in the total operating costs. 

[0021] The rate of ageing e{x) is computed from models associated with various plant 
components, such as the following models, which are well know to one versed in the art: 
[0022] Microcrack growth rate models for components operating at elevated 
temperatures, such as boiler pressure parts, steam pipes and steam turbine rotors. See for 
example "Life Assessment Technology for Fossil Power Plants", R. Viswanathan, Proc. 
American Power Conf., Vol. 57-111: 1749-1768, 1995 and "A Lifetime Prediction Model 
for Thermomechanical Fatigue based on Microcrack Propagation", M.P. Miller et. al, 
ASTM STP 1186, 1993, pp. 35-49. 

[0023] Corrosion, oxidation, and erosion models based on extensions to microcrack 
growth models as described in "Corrosion Fatigue: Mechanics, Metallurgy, 
Elctrochemistry, and Engineering", T. Crooker and B. Leis editors, ASTM 1983. 
[0024] Mechanical and thermomechanic stress models to determine stresses created in 
plant components, e.g. stresses in a steam turbine rotor shaft due to load oscillations and/or 
steam temperature changes. 



[0025] The rate of ageing e(x) represents the plant lifetime consumption or an asset value 
degradation expressed as a financial value under operating conditions caused by the control 
values u{x) . The entries in e(x) correspond to the aging rate for each plant component and 
are, for example, computed directly from 

R da 

a crit dr ' 

where a crit is a predefined critical crack length for said component, a is the crack length 

computed using the models described above and dependent on U(x) , and, R, is the 

component, cost, repair, or maintenance cost. 

[0026] In the inventive method, the following steps are performed 

[0027] a) The optimization module 1 determines the trajectory of plant control values 

u(x) from the present time t for a predetermined duration T and transmits said plant control 

values U(x) to the simulation module 2. The first time that this step is executed, plant 

control values are e.g. set to constant values for the duration Tor are determined from 

stored standard trajectories corresponding to a given operating condition such as a plant 

startup or a load change. Later executions of this step preferably use a variation of the 

plant control values u(x) of an earlier step. 

[0028] b) The plant model incorporated in the simulation module 2 determines from 
the trajectory of control values u(x) , by simulation of the power plant behavior, a 
trajectory of the state of the power plant. From said state, corresponding trajectories of 
costs c(x) for operating the plant and trajectories of revenues q(x) from selling produced 
output are determined, as well as the rate of ageing e(x) of plant components. In a 



preferred embodiment of the invention, trajectories of plant output values p(x) are 
determined as well. 

[0029] c) A total operating cost JfuJ is computed as an integral of an objective 
function that comprises the rate of ageing e(x) and a difference between costs ?(t) and 
revenues q(x) . In a preferred embodiment of the invention, this total operating cost JfuJ : 
defined as the functional 



where R l5 R 3 are weighting matrices. The weighting matrices determine a relative 
influence of the magnitude of each component of each of the different vectors 
q(x),c(x),e(x). 

[0030] In a preferred embodiment of the invention, the above objective function also 
includes a difference between plant output values p(x) and demand values d(x) . In this 
embodiment, the total operating cost JfuJ is preferably defined as 



where R l5 R 2 , R 3 are weighting matrices. The terms in the objective function are 
interpreted as follows. The "d(t) -p(vy term represents a penalty associated with not 
satisfying the required power output. The u c(v) - q(r)" term represents costs minus 
revenue, i.e., an economic performance indicator. 




J[u] = f R x e(x) + R 2 {d(x)-p{x)f + * 3 [<?(t)-?(t)] dx, 
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[0031] d) Steps a) through c) are repeated iteratively, with the optimization module 
varying the trajectory of plant control values u:[t, t + T]~R" until an optimized lower total 
operating cost J[u] is arrived at. As the functions q(x),q(x), c(x), e(x) depend upon plant 
control values u(x) , the optimization routine finds a strategy u:[t, t + T]»R n that minimizes 
the total operating cost JfuJ. 





[0032] 


Examples for units are, 




[0033] 


for the rate of ageing e : $ per second 


"y 


[0034] 


for costs c or revenue q : $ per kW or $ per kg steam per second 


[0035] 


for Rl : dimensionless 




[0036] 


for R2: $/kW A 2 or $(s/kg) A 2. 




[0037] 


for R3: kW or kg per second 




[0038] 


for demand d or production p : kW or kg per second 




[0039] 


In another preferred embodiment of the invention, the objective function is 



J[u] = f R,e(x) + R 2 [d(x) -p(x)f + R 3 c(x) -R 4 q(x) dx, 



where costs and revenues are weighted separately, using an additional weighting matrix R 4 , 
which allows for the cost and revenue vectors c(x) , q(x) to have different dimensions. 
[0040] The optimization step c) involves the problem of how to find the optimal solution 
efficiently, since the search space determined by the plurality of plant control values and 
the duration T is large. Standard optimization procedures may be used, or dedicated 
solutions as described, for example, in US 5,347,466, US 5,886,895 or in "The 
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Generalized Unit Commitment Problem", Ross Baldick, IEEE Transactions on Power 
Systems, 10(l):465-475, February 1995. 

[0041] In a preferred embodiment of the invention, at time t, the control strategy, i.e. the 
future values u:[t, t + T]*R n corresponding to optimal plant operation, is displayed 
graphically to an operator, providing him with guidance for controlling the plant. 
[0042] Figure 2 shows a block diagram corresponding to another preferred embodiment 
of the invention, in which, at time t, the control value u(t) is applied to the actual power 
plant 3, resulting in automatic optimized plant operation. Preferably, measured plant 
values y(t) of the plant are input to the simulation module 2 in order to update the 
simulation model. Appropriate modeling and model update techniques are well known 
from the field optimal control. In particular, receding horizon control techniques might 
result useful. The power plant is, for example, a gas turbine, a combined cycle power 
plant, or of a boiler-turbine-generator type. 

[0043] A computer program product according to the invention comprises a computer 
readable medium, having thereon: computer program code means to make, when said 
program is loaded in a computer that is operationally connected to a power plant, the 
computer execute the method according to one of claims 1 to 6. A computer program 
according to the invention comprises said computer program code. 
[0044] A power plant optimizing system that optimizes plant control values according to 
the invention comprises 

[0045] a) an optimization module 1 that is configured to determine a plurality of 
trajectories of plant control values u(x) and to determine a total operating cost J[u] as an 
integral of an objective function that comprises a difference between costs c(t) and 
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revenues q{x) that are associated with operating the power plant according to the plant 
control values u(x) , and 

[0046] b) a simulation module 2 that is configured to determine, from the trajectory of 
plant control values u(x) , by simulation of the power plant's behavior, corresponding 
trajectories of costs c(x) for operating the plant and trajectories of revenues q(x) from 
selling produced output. 

[0047] The simulation module 2 is configured to determine, from the trajectory of plant 
control values u(x) , a trajectory representing a rate of ageing e(x) of plant components. 
The objective function comprises said rate of ageing e(x) . 

[0048] In a preferred embodiment of the power plant optimizing system, the simulation 
module 2 is configured to determine, from the trajectory of plant control values u(x) , plant 
output values p(x) , and the objective function comprises a difference between plant output 
values p(x) and associated demand values d(x) . 

[0049] The invention allows the operator or an automatic controller to incorporate not 
only short term financial considerations into a plant operation strategy, but also long term 
financial considerations that involve maintenance and replacement costs and their 
dependency on operating conditions. 
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List of designations 




[0050] 


1 


optimization module 


[0051] 


2 


simulation module 


[0052] 


3 


actual power plant 


[0053] 


u(v) 


control values 


[0054] 




production 


[0055] 


C(T) 


costs 


[0056] 


d{z) 


demand 


[0057] 




revenue 


[0058] 




rate of ageing 


[0059] 




measured plant values 


[0060] 


R l5 R-25 R35 


weight matrices 



